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Abstract 

We relate the fragmentation functions of partons into mesons and baryons in the 
framework of recombination of shower partons. The results are in reasonable agree- 
ment with the data. The implication is that the meson and baryon fragmentation 
functions are not independent when hadronization of the shower partons are taken 
into account. The conclusion therefore closes a conceptual gap in the system of frag- 
mentation functions whose Q 2 evolution has been more extensively studied than their 
interrelationship. 

PACS numbers: 25.75.Dw 



1 Introduction 

Hadron production at intermediate pt in heavy-ion collisions (HIC) has been successfully 
described in the recombination model (RM) [1, 2, 3]. In Ref. [1] the hadronization process is 
shown to be dominated by the recombination of thermal and shower partons when pt is in 
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the range 3 < pr < 8 GeV/c. The distributions of those shower partons are derived by fitting 
the fragmentation functions (FFs) in the RM [4]. The shower partons considered are not to 
be identified with those due to gluon radiation at very high virtuality that are calculable in 
pQCD, since the ones we consider are defined by their recombination within a jet to form a 
hadron at the hadronization scale, and are not calculable in pQCD. In this paper we revisit 
the problem of fragmentation in the framework of recombination, and calculate the FFs for 
proton and A using the shower parton distributions (SPDs) that have been determined in 
Ref. [4] from fitting the FFs for mesons. 

The determination of FF for proton in a jet that can agree with data would be a significant 
achievement, since it is the first time a FF is obtained by calculation, instead of by fitting. 
It should be stated from the outset that we are not concerned with Q 2 evolution here, since 
that part of the problem is well handled by pQCD [5]- [8]. All evolutions have to start from 
somewhere, and the starting distributions have to be determined by fitting the data. In Refs. 
[6, 7] the FFs have been parameterized for a variety of initiating partons and final hadrons, 
and have all been determined independently apart from trivial identities. To our knowledge 
no one has attempted to relate one FF to another. Our attempt to relate D 7r (z) and D' p (z), 
FFs for 7r and p, respectively, is feasible because our emphasis is on the hadronization of 
the shower partons. Thus a successful connection among the FFs is an affirmation of the 
soundness of the recombination formalism that we employ, and establishes an important link 
among the multitude of FFs. 
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2 Gluon Fragmentation into Proton 

We first recall how we have treated the fragmentation of gluon into pion in the framework 
of recombination. The basic equation is [4] 

xD* G {x) = I ^^{G(x 1 ),G(x 2 )}R 7r (x 1 ,x 2 ,x) , (1) 

J X\ X 2 

where Dq(x) is the FF of gluon into pion, and R n (xi,x 2 ,x) is the recombination function 
(RF) for qq — > n that is given in Ref. [1, 9]. G(x,j) is the SPD of q or q in a gluon jet, and 
the symmetrization symbol in Eq. (1) means 



{G(x 1 ),G(x 2 )} = - 



G( Xl )G ( j-^-) + G{x 2 )G (j Xi 



(2) 



XxJ \l-x 2 

Since Dq(x) is known from the parametrization by BKK [5] determined by fitting various 
data, we can solve for G(z) in Eq. (1), not analytically, but by fitting the FF at a chosen 
Q 2 . The form for G(z) used is 

G(z) = Az a (l - z)\l + cz d ) . (3) 

For Q 2 = 100 GeV 2 the values of the parameters are A = 0.811, a = -0.056, b = 2.547, 
c = —0.176, and d = 1.2 [4]. The evolution of these parameters for other values of Q 2 has 
been considered elsewhere [10]. It should be emphasized that G(z) is defined by Eq. (1) in 
the context of parton recombination, and is the distribution at the hadronization scale of a 
light quark (or antiquark) in a gluon-initiated jet. The q and q at x\ and x 2 are to recombine 
to form a pion at x, which is constrained by R n (xi, x 2 , x) to be x = X\ + x 2 . 

We now advance the point of view that if Eq. (1) makes sense in describing FF in terms 
of the recombination of shower partons, then the same formalism should make possible a 



calculation of the FF of gluon into proton. Indeed, we need only generalize Eqs. (1) and (2) 
in the self-evident way 



xDg(x) = I — — — {G(x 1 ),G(x 2 ),G(x 3 )} R p (x 1 ,x 2 ,x 3 ,x) , (4) 
J X\ x 2 x 3 



where 



{G(x 1 ),G(x 2 ), G(x 3 )} — ^ 



(5) 



there being six terms in the permutation of x 1: x 2 and x 3 . Dq(x) is an initial form of the 
FF of gluon into proton, on which modification will be made in the following. R p is the RF 
for the proton [1] 



/ x u x 2 \ a+1 (x 3 y +1 (Xx +x 2 + x 3 x 
R p (x 1 ,x 2 ,x 3 ,x) = g at g p [—j-) {-) 5 1^ 1) , (6) 



where g st = 1/6 and g p = [B(a + 1, a + (3 + 2)B(a + 1, (3+ 1)] 1 , B(a, (3) being the Euler-beta 
function. The parameters a and (5 are [11] 

a = 1.75, /5 = 1.05 . (7) 

We do not regard xDq(x) as the final form for the FF because we have not so far 
taken into account the momentum carried by the antiproton that must accompany the 
production of a proton in a gluon jet. To treat the pair production of pp, we have to 
consider the correlation between the dibaryon produced. If D P q(x) denotes the pp FF, then 
upon integration over x 2 we should have on general grounds 

dx 2 D%( Xl ,x 2 ) = D p G ( Xl ) [ dx 2 D%{x 2 ) (8) 
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which does not imply factorizability of x 2 ). We shall apply the momentum constraint 

x\ + x 2 < 1 as the only correlation between p and p, and write 



D™( Xl ,x 2 ) = X - 



(9) 



Commensurate to Eq. (9), we approximate D G (x 2 ) in Eq. (8) by D'§(x 2 ). Then D G (x±) can 
be evaluated from 



jdx 2 {D'S(x 1 ),D'j (x 2 )} 
J dx 2 D%(x 2 ) 



dim = - — L fJ _ V,_r J • (10) 



The result for xD p G (x) is shown in Fig. 1 by the solid line. Also shown in that figure in 
dashed line is the FF obtained by KKP [7] by fitting the data in NLO. The agreement over 
5 orders of magnitude is very good, considering the fact that no free parameters have been 
used in the calculation. 



3 Other Fragmentation Functions 

In this section we consider two other FFs into baryons, first g — > A, and then u — > p. For 
the former Eq. (4) is replaced by 

xD'i(x) = [({[—) {G(x 1 ),G(x 2 ),G s (x 3 )}R A (x 1 ,x 2 ,x 3 ,x) , (11) 

J \i=l x i / 

where G s (xs) is the SPD of a strange quark in a gluon jet, and R\(xi, x 2 , £3, x) is the RF 
for A. In Ref. [4], G s {z) is parameterized as in Eq. (3), but with A = 0.069, a = —0.425, 
b = 2.489, c = —0.5 and d = 1.1. _Ra is given in Ref. [12]; it has the same form as Eq. (6) 
but with g st = 1/4, = g p , a = 1 and (3 = 2. Thus D'q(x) can be calculated using Eq. (11) 
and then modified as before to give D G (x) by taking the momentum of A into account, as 
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done in Eq. (10). The result is shown by the solid line in Fig. 2. For comparison, we show 
in the same figure the LO and NLO distributions of g -> A at Q 2 = 100 GeV 2 given in Ref. 
[13]. Although our result is a little lower in the mid-x region, the agreement is nevertheless 
acceptable; in fact, it agrees better with the LO result. 

Next, we consider proton production in a jet initiated by a u quark. In that case valence 
and sea quark distributions must be distinguished. In Ref. [4] we have used K NS and L to 
denote the valence and sea quark SPDs, respectively; they are parameterized as in Eq. (3) 
with the parameters for K NS (L) : A = 0.333 (1.881), a = 0.45 (0.133), b = 2.1 (3.384), c = 
5.0 (—0.991), and d = 0.5 (0.31). The recombination formula is now 



xD*(x) =[({[—) {K(x 1 ),L(x 2 ),L(x 3 )}R p (x 1 ,x 2 ,x 3 ,x) , (12) 
J \i=i x 't I 



where 



K( Xl ) = K NS (x 1 ) + L( Xl ) . (13) 

For p production in the same jet valence quark does not contribute, so for D'^{x) the sea 
quark distribution L(xi) replaces K{x\) in Eq. (12). The resultant D'£ and D'£ are then 
substituted into 

DP(r , _ Sdx 2 {D^x 1 ) i m(x 2 )} 
u[l) ~ jdx 2 D'S{x 2 ) ■ (i4j 

The FF thus obtained is shown by the solid line in Fig. 3. The parametrization of the data 
by KKP [7], shown in the dashed line, is, however, unreliable because an invalid assumption, 
DP(z) = 2D p d (x) ) has been used in the analysis. While we await an improved analysis of the 
data, the solid line in Fig. 3 may be regarded as our prediction of the FF for u quark into p. 
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4 Conclusion 

We have calculated the FFs into baryons, using the SPDs determined by fitting the FFs 
into mesons. The results are in good agreement with data for g — > p and g — > A. For 
■u — > p there is no reliable analysis of the data, so our result is a prediction. To be able to 
relate meson and baryon FFs is an attribute of our formalism that treats fragmentation as 
a recombination process. Since the shower partons are defined at the hadronization scale, 
their distributions cannot be calculated from first principles. Nevertheless, they have been 
determined by fitting FFs for mesons, inasmuch as the FFs themselves have been determined 
by fitting the data. What we have achieved in this paper is to elevate the significance of 
the formalism that connects previously unrelated FFs. That not only confirms the internal 
consistency of the formalism, but also renders additional support to the reliability of the 
procedure employed earlier in the study of the contribution of thermal-shower recombination 
to hadron production (mesons and baryons) at intermediate pt in HIC. 

The experience gained through this work on baryon production builds the basis for future 
investigation of correlations that involve baryons produced at intermediate px in HIC. Data 
on such correlations from experiments at RHIC already exists [14]. Whether a baryon is 
the trigger particle or an associate particle, it is clear that the recombination model is well 
positioned to address the problem of hadronization of various partons that arise from different 
sources (hard and soft) in HIC. 



7 



Acknowledgment 

We are grateful to W. Vogelsang for helpful communication on fragmentation functions. 
This work was supported in part, by the U. S. Department of Energy under Grant No. 
DE-FG02-96ER40972 and by the Ministry of Education of China under Grant No. 03113. 



8 



References 

[1] R. C. Hwa and C. B. Yang, Phys. Rev. C 67, 034902 (2003); 70, 024905 (2004). 

[2] V. Greco, C. M. Ko, and P. Levai, Phys. Rev. Lett. 90, 202302 (2003); Phys. Rev. C 
68, 034904 (2003). 

[3] R. J. Fries, B. Muller, C. Nonaka and S. A. Bass, Phys. Rev. Lett. 90, 202303 (2003); 
Phys. Rev. C 68, 044902 (2003). 

[4] R. C. Hwa, and C. B. Yang, Phys. Rev. C 70, 024904 (2004). 

[5] J. Binnewies, B. A. Kniehl and G. Kramer, Z. Phys. C 65, 471 (1995). 

[6] B. A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B582, 514 (2000). 

[7] B. A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B597, 337 (2001). 

[8] S. Kretzer, Phys. Rev. D 62, 054001 (2000). 

[9] R. C. Hwa and C. B. Yang, Phys. Rev. C 66, 025205 (2002). 
[10] Z. Tan and C. B. Yang, nucl-th/0502008. 
[11] R. C. Hwa and C. B. Yang, Phys. Rev. C 66, 025204 (2002). 
[12] R. C. Hwa and C. B. Yang, OITS-773 (to be submitted for publication). 
[13] D. de Florian, M. Stratmaan and W. Vogelsang, Phys. Rev. D57, 5811 (1998). 
[14] S. S. Adler et al, (PHENIX Collaboration), Phys. Rev. C 71, 051902 (2005). 

9 



Figure Captions 

Fig. 1. Fragmentation function for gluon into proton. Solid line is our calculated result; 
dashed line represents the fit of the data in NLO [7]. 

Fig. 2. Fragmentation function for gluon into A. Solid line is our calculated result. The line 
with circles represents the fit of the data in NLO, while the dashed line with squares 
is in LO [13]. 

Fig. 3. Fragmentation function for u quark into proton. Solid line is our calculated result; 
dashed line represents a questionable fit of the data in NLO [7] (see text). 
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